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Abstract. The expression for polarized electric dipole moment of well-deformed reflection asymmetric 
nuclei is obtained in the framework of liquid-drop model in the case of geometrically similar proton and 
neutron surfaces. The expression for polarized electric dipole moment consists of the first and second 
orders terms. It is shown that the second-order correction terms of the polarized electric dipole moment 
are important for well-deformed nuclei. 

PACS. 23.20.-g Electromagnetic transitions - 23.20.Lv Gamma-transitions and level energies 



43 ! 1 Introduction 



Reflection asymmetric deformation of nucleus induces the 
proton- neutron redistribution. As a result, the proton or 
neutron density distributions became slightly polarized 
and reflection asymmetric in the nuclear volume. Due to 
such density polarization the position of proton center of 
mass is shifted relatively the nuclear center of mass; there- 
fore reflection asymmetric nuclei have the polarized elec- 
tric dipole moment (PEDM). 

The PEDM of nuclei with quadrupole and octupole 
surface deformations was firstly obtained by V. M. Struti- 
nsky in 1956 [Tj in the framework of liquid-drop model. A 
short time later A. Bohr and B. R. Mottelson was evalu- 
ated the PEDM in the same model [2J, but Strutinsky's 
derivation is the correct one [3]. The PEDM was found 
for non-axial nuclei with quadrupole and octupole defor- 
mations in Refs. [U[S]. Note the PEDM discussed in Refs. 
[1,2,3,4,5 is only related to the proton- neutron polariza- 
tion in the volume of nuclei with quadrupole and octupole 
surface deformations. 

However the proton-neutron density polarization in 
the nuclear volume induces the variation of proton and 
neutron radii and, therefore, leads to the corresponding 
surface contribution into the PEDM. The expression for 
PEDM with volume and surface contributions was derived 
in Ref. [5] in the framework of the droplet model for axial 
nuclei with the proton radius R p (9) in the form 



R P (0) 



R 



Op 



= F(9) = 



(1) 



Here Rq p is the proton radius of spherical nucleus, j5g is 
the deformation parameter and Yg {9) is the spherical har- 
monic function. Both spherical and deformed nuclei have 



a neutron skin of constant thickness in the framework of 
the droplet model, see Ref. [B] and papers cited therein, 
therefore the neutron radius R n {9) of deformed nuclei 
is not proportional the proton one R p (9) in the droplet 
model. Due to this expression for the PEDM obtained 
in the droplet model consists of the volume and surface 
charge redistribution contributions as well as the contri- 
bution related to the neutron skin thickness [BJ[7]. The 
neutron skin thickness contribution arises precisely from 
the non-coincidence of the centers of mass of a uniform 
skin and of the volume it encloses [6l[7]. 

The expression for PEDM for nuclei with geometrically 
similar proton and neutron surfaces, i.e. when proton and 
neutron radii have the same angular dependence 



R P {9) _ Rn{0) 



R 



Op 



Ror, 



F(0), 



(2) 



was obtained in Ref. [5] . Here i?on is the neutron radius of 
spherical nucleus. The neutron skin thickness depends on 
9, when radii of the proton and neutron surfaces are pro- 
portionally to each other. The PEDM consists of the vol- 
ume and surface charge redistribution contributions only 
in this case [8], because the neutron skin thickness con- 
tribution equals zero for geometrically similar proton and 
neutron surfaces, see for details [TUSH] and Sec. II. The 
neutron skin center of mass coincides with the nucleus the 
ones [5J|5]. 

The expressions for volume and surface parts of the 
PEDM for non-axial nuclei with arbitrary multipole defor- 
mations and geometrically similar proton and neutron sur- 
faces are given in Ref. |10) . The equilibrium shapes of some 
nuclei are non-axial reflection asymmetric |11) . Moreover 
PEDM can be arisen at non-axial reflection asymmetric 
surface vibrations |10) . 
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We emphasize that the PEDM obtained in the first 
non-zero order on multipole deformations of nuclear sur- 
face is proportional to /3f/3^+i and all expressions for the 
PEDM presented in Refs. [TPPP1|31|51|7||5P^ are de- 
rived in this approximation. 

Numerical study of the PEDM in well-deformed nuclei 
in Ref. [12] shows that the first approximation for PEDM 
is strongly underestimated the numerical one. Moreover 
the difference between the numerical and first-order values 
of PEDM increases with values of deformation parameters 
strongly [T2"] . 

The values of PEDM have been also evaluated in the 
frameworks of various semi-microscopic or microscopic ap- 
proaches, see Refs. [9l fT2"I[T3T 1 4 II 1 5 [HB] and papers cited 
therein. 

The nuclei with quadrupole and octupole deforma- 
tions, El transitions and the PEDM are studied inten- 
sively recently p^rT^rT^f20lf2Tlf2^f23lf24lf25lf26| . The PEDM 
plays important role in various phenomena of well-deformed 
reflection asymmetric nuclei. Thus Karpeshin has shown 
that well-deformed fission fragments of such shapes formed 
during prompt fission give rise to both the anomalous El 
internal conversion |27| and the prompt gamma radiation 
[251125] related to the PEDM. The left-right asymmetry of 
fission induced by polarized neutrons can be also linked 
to the PEDM [30]. The El transitions possibly linked to 
octupole vibrations around super-deformed shape can be 
also enhanced by the PEDM [3T1I521I531IM1 I33] . Strong El 
transitions related to the low-energy shape oscillations of 
negative parity in the first and second (isomeric) minima 
in actinides are also connected to the PEDM [36 . 

However application of expression for the PEDM ob- 
tained in the first order for well-deformed nuclei is ques- 
tionable as pointed by Skalski [12] . Therefore it is desir- 
able to obtain expression for the PEDM in the next order 
at least, which is the second order approximation for the 
PEDM contained terms proportional to Such ex- 

pression should be helpful and practical for description of 
various nature El transition in well-deformed nuclei. 

The PEDM evaluated in the framework semi-microscopic 
approaches [13 , 14 consists of macroscopic and microscopic 
contributions. The microscopic contribution of PEDM is 
evaluated without applying the perturbation approach on 
the surface deformation parameters, while macroscopic 
the one is evaluated by using expressions obtained in the 
first non-zero order on multipole deformations. Therefore 
more accurate expression for the macroscopic part of the 
PEDM improves the accuracy of the PEDM evaluated in 
the framework of semi-microscopic models. 

It is well-known that shell effects arc reduced in heated 
nuclei. Due to this expression for PEDM has reliable ac- 
curacy for highly excited fission fragments. The shape of 
fission fragments after rupture has appreciable reflection 
asymmetry. Microscopic calculation of PEDM in heated 
nuclei has not done up to now. Therefore expression for 
the PEDM obtained in the next order at least can be useful 
for evaluation of various effects related to dipole moment 
of fission fragments. 



Various proposals on neutron skin are discussed re- 
cently. The neutron skin of permanent thickness in de- 
formed nuclei is widely applied in the framework of the 
droplet model [SE] and other papers related to this model. 
The neutron skin with thickness depended on 9 and re- 
lated to well-known and widely-used relation for the pro- 
ton and neutron radii (2) is also very common in nuclear 
physics; see, for example, the description of the proton and 
neutron mean field radii and etc. HESMHlEIlEHlIMlIinillll 
42 , 43 , 44 . We emphasize that the shapes of potential and 
density distributions are geometrically similar due to a 
consistency between the density distribution and the cor- 
responding mean-field potential 9 . More complex approx- 
imation to the neutron skin thickness related to different 
neutron versus proton deformations of nuclei is discussed 
in Ref. 45 , 46 . The thickness of neutron skin has complex 
angular dependence on 9 in such case. 

Diverse approximations for neutron skin shape are rea- 
sonable for small and medium deformed nuclei close to 
the beta-stability line. The ratio (2) between the proton 
and neutron radii is, probably, discussable for extremely- 
deformed and/or neutron- rich nuclei as pointed in Ref. [7]. 
Current experimental studies of neutron distribution on 
the surface of nuclei devote to the neutron skin thickness 
in spherical nuclei mainly, see Refs. [4"Tll4"8] and papers 
cited therein. Available experimental data [47,48 cannot 
support firmly any of these approximations on the neutron 
skin thickness in deformed nuclei. 

Relationship between proton and neutron surfaces or 
potential described by Eq. (2) is very widely used in nu- 
clear physics for small, medium and even well-deformed 
nuclei [9lfT^fT^[3Tlf3^[3^[40l[4Tl [4 ^ [43 l [44 ] . Therefore, in Sec 
II we derive expressions for the PEDM for the geomet- 
rically similar proton and neutron surfaces, which takes 
into account the first (oc PtPt+i) and second (oc fitPl'fit") 
orders contributions. Discussion of obtained expressions, 
numerical results and conclusion are given in Sec. III. 

2 Model and expression for PEDM 

Let us consider the axial nucleus with proton and neutron 
radii described by Eq. (2). There are no any density polar- 
izations in spherical nuclei, therefore the equilibrium neu- 
tron and proton density distributions in deformed nucleus 
can be presented as p n — po n + Sp n and p p = po p + Sp p . 
Here p 0n = 37V/(47ri?Q„) and p 0p = 3Z/(4ttR% p ) are the 
equilibrium neutron and proton densities in spherical nu- 
cleus, 5p n and 5p p are the variations of neutron and proton 
densities induced by surface deformation, Z and N are the 
numbers of protons and neutrons in the nucleus. 

Due to high value of the nuclear matter incompress- 
ibility the total nuclear density p = p n + p p in the nuclear 
volume is constant p = po n + po p , therefore Sp n = —Sp p , 
see also [UH]. 

We should take into account that the numbers of pro- 
tons and neutrons in deformed nucleus are, respectively, 
Z and N; and the center of mass must lie in the plane of 
mirror symmetry of the nucleus [T][H]IS] , because the reflec- 
tion asymmetric nuclear shapes are coupled by sub-barrier 
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tunnel transition. These two conditions can be easy ful- 
filled by introduction of auxiliary monopole /3q and dipole 
Pi deformations, i.e. 



Rp(8) _ R n (0) 
Ron 



R, 



Op 



F{6) + P o Y oo (0) + PiY 10 (d) 



L 



= f(9) = 1 + > JtYtaiS). (3) 



1=0 



The values of Pq and Pi are, correspondingly, determined 
by equations 

JdV P f = JdV^ = l£de sin {9)f{9f = 1, (4) 

3 

dV r cos (9) {p p + p 0n ) = o ( ZR 0p + NR 0n ) 



x / d6 sin (9) cos (9) f(6f = 0.(5) 
Jo 

For the sake of simplicity we take into account the 
most important multipole deformations of nuclear surface 
Pi, p3, Pa, /?5, Pq. The expressions for Po and Pi taken into 
account all quadratic and cubic terms on p2, P3, Pi, P5, Pe 
can be directly obtained from Eq. (4) and (5), however cor- 
responding equations are cumbersome and therefore not 
presented here. 

The PEDM is defined as 



(0) 



D = e dV r cos (9) pp. 



Due to deviation of the nuclear surface from spheri- 
cal form there are variation of the proton density into the 
nuclear volume Sp p (r). The variation of nucleon density 
in nuclear volume induces the deviation of the proton ra- 
dius 5Rp(9) from the equilibrium position on the nuclear 
surface. The proton radius variation induces the proton 
density variations in the volume SR P (9)AS, where AS is 
the element of surface square. Therefore the PEDM in re- 
flection asymmetric nuclei with axial symmetry is related 
to the redistribution of protons relatively neutrons into 
the nuclear volume and on the nuclear surface, see also 
0, 



D = D V +D S , 



(7) 



where 



D v e J dVr cos (6) [po p + Sp p ] — e J dVr cos (9)Sp p 

r-* pRopf(e) 
= 2ire d9 sin (0) cos (9) / dr r 3 Sp p , (8) 

Jo Jo 

D s « e J dS R p {9) cos (9) p 0p SR P {9) 

= — - / dd sin (9) cos (0) 
2 Jo 



3Ze r 
A) 

xf(9)6Rp(9). 



1+ rm 

\f(8) 



1/2 



(9) 



Here Eq.(5) is taken into account at simplification of Eq. 
(8) and f(6) = So, the volume part of PEDM is 

related to the volume integral and density variation in the 
nuclear volume 5p p while the surface part of the one is 
determined by the surface integral and the proton radius 
variation SR P (9). 

The proton (or neutron) density variation induced by 
surface deformation produces additional pressure on the 
free nuclear surface. Due to this pressure the position of 
corresponding surface is slightly shifted. Both the surface 
symmetry energy and Coulomb force counteract the sur- 
face shift and neutralize the additional pressure on the 
free nuclear surface induced by density variations, see for 
details 6 , 8 :,[TD] . Normal to the surface variation of the pro- 
ton radius is defined by the boundary condition [61I81IT0]. 
which equalizes the normal to surface pressures induced by 
density fluctuations, neutron-skin stiffness and Coulomb 
interaction, and equals to 



SR P (9) 



_N 3ei? 
~ ~A 8QA 1 / 3 

N 3ei? 



<f>{R p {9)) 



J dS<t>(Rp(9)) 



A 8QA 1 / 3 



y(Rp(9)) - 



JdS 

fdS<p(R p (9)) 
JdS 



(10) 



where Q is the neutron-skin stiffness coefficient [MTU] , <fi = 
ip—Tp, (p(r) is the Coulomb potential related to the protons, 
_ j dVip 
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dV 



is the average potential value in the nucleus 



and A = Z + N. Note that Z5R P {9) + N5R n {9) = 0, 
because the center of the mass must lie in the plane of 
mirror symmetry of the nucleus, i.e. 



dV r cos (9) p(r) 



= J dV r cos (9) [p p + Sp p + p 0n + Spn] 
+ J dS cos {9) [Rp(9)popSR p + R n (9)p Qn SR n ] 
= J dV r cos (9) [Sp p + Sp n ] 
+ J dS cos (9) [Rp(9)popSR p + R n (9)p Qn SR n ] 
dS cos{9) [Rp{9)p Qp 5Rp + R n (9)p 0n 5R n ] 
dS cos (9) Rp(9)[pop8Rp + p 0n 8R n ] = 0. 



Here Eq. (5) and condition Sp p — —5p n are taken into 
account. 

If we know Sp p (r) and <p(r) than we can evaluate the 
PEDM using Eqs. (7)-(10). Let us find Sp p (r) and ip(r) in 
the framework of liquid-drop model. The energy density 
functional, which is described density distribution in the 
nuclear volume, can be written in a simple form [T1I51I10) 



t w — a v p + J 1- eppip 

P 

(p-2ppf 
= -a v p + J h eppip, 



(11) 
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where — a v is the bulk energy per nucleon in symmetric 
nuclear matter and J is the volume symmetry energy. Note 
that the energy density functional of the droplet model 
contains the dilatation term [BJ, but parameter L related 
to the dilatation term is equal zero in recent parameter 
set of the droplet model [B] , therefore we neglect dilatation 
term here. The energy of nucleus E is related to the energy 
density functional E = J dV£. The equation determined 
the equilibrium distribution of the charge into the nuclear 
volume can be obtained by variation of the energy 

SE = 6 j dV [£ - Xpp] (12) 

= dV [-4J0 - 2p p )/p + eip~ (a v + A)] Sp p 



on 5p p with the additional condition conserved the number 
of protons in the nucleus. As the result, we get 



8Jp p = -p(etp - 47 - A'), 



(13) 



where A' = a v + A and A is the Lagrangian coefficient 
related to the additional condition. The solution of this 
equation is 



Pop = P o + 



1 



A' 



Sp p 



2 8 J 8 J ) 

, (\ A' eTp\ 
(po P + Pon) + 

-ep(ip - Tp) _ ~e(po p + po n )((p - Tp) 



8J 

— 3eA(ip — Tp) 

327Ti?3 p J 



8.7 



-3eA(p 
32nR 3 p J' 



(14) 



(15) 



Note that / dVdpp = 32n 3 ^ j U dV ^ ~ S 5 )] = °- We ne " 
glect by the difference between Ro p and Ro n at evaluation 
of volume quantities. 

The electric potential can be found by using the Pois- 
son equation 



V cp = Airepp 



(16) 



Substituting (15) into (8) and taking into account (5) 
we get 



-3e 2 A 
16 JRlp Jo 



rRopf(e) 

d6 sin (0) cos (0) / drr 3 ^ 
Jq 



-3e 2 A 



WJR 3 0p jo 



rRo P f(e) 

d0 sin (0) cos (0) / drr 3 ip(r). (17) 

Jo 



Using (10) and approximation sa 4 ; se e Ref. [BJ, we 
rewrite (9) into the form 



9Ae 2 R 



d0g{9) cos(0)f(0) 

£ dff g(0'MR p (0')) 
Iod0'g(0') 



(18) 



1/2 



So, the 



where g{0) = sm(9)f 2 (0) [l + (f'(d)/f(0)) J 

volume (17) and surface (18) parts of PEDM are deter- 
mined by the Coulomb potential ip(r). 

The Coulomb potential of deformed nucleus is 



ip(r) 



dV 



PpW e f dV !*z+M*l . (19) 



This potential satisfies to Eq. (16). 

It is possible to find potential tp(r) by applying the per- 
turbation theory to Eqs. (14), (15) and (19). We expand 
the potential and the variation of proton density into the 
perturbation series 

<p(r)= V >°(r) + <p 1 (r) + l p 2 (r) + ... ) (20) 
Sp p (r) = 6p° p (r) + 6 Pp (r) + Sp 2 p (r) + ... , (21) 

where the superscript corresponds to the number of per- 
turbation approach. Corresponding solution for the La- 
grangian coefficient is variation of proton density into the 
perturbation series 



A = -a v + 8J[Z/A - 1/2 + e/(8J)(<p°(r) 
+</(r) + <^ 2 (r) + ...)]. 

Substituting perturbation series (20)-(21) into Eqs. (15) 
and (19) we get 



-ep4> k {r) -ZeA<t> k {r) 



8.7 



32ttR 3 p J 



, for k > 0, (22) 



/(r) 



dV 



5p k p -i(r') 



for k > 1. 



(23) 



and 



dV- 



Po P 



eZ 
Rq p 



dV 



4^ p |r-r' 



eZ \ - 

Ro~p^ (2£+l) J 



d0' $m(0')Ylo{0') 



dr'- 



(r') 



>\t+2 



Ra P f(e') 



dr'- 



{r'Y- l {Ro P f 



(24) 



where 4> k (r) = ip k (r) — ip k (r). 

Using Eqs. (22) and (23) we get the recurrent equation 
for <p k (r) at k > 1 



/(r) 



-3e 2 A 



32irR 3 p J 



dV 



fc - 1 (r / ) 



ElS^ / d0>sm(O')Y; o (0>) 



R Qp J Y 4(2* +1) Jo 



dr 



,(r')^ + V fc_1 (r') 



ri+l 



(Ro P y 



dr 



(r'Y-HRop)' 



(25) 
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which determines the potential with any necessary degree 
of accuracy. As a result, we can evaluate the volume and 
surface contributions ofPEDM using Eqs. (17), (18), (20), 
(24) and (25). 

The macroscopic PEDM can be written as 



A, 



where 



Ai + D v 



Ai = 



A, 



i 3 AZ 



Ai + D s 



9/3 2 /3 3 , H/3 3 /34 



A 



.(26) 



5G 



41/3 4 /3 5 



D 



264 

15e 3 4 2 / 3 Z 



si 



35 105a/7 

441/3 5 /3 6 " 
11 715a/143 

9/3 2 /3 3 , ll/3 3 /3 4 



SttQ |_56a/35 105a/7 

41/3 4 /3 5 , 441/3 5 /3 6 ' 



are the volume and surface hrst-order contributions, 



(27) 



(28) 



A 



l AZ 



!U ~~ 7T 3 / 2 J 

65685/3 2 fe^ 6 
44044 a/455 



3/3 2 /3 3 , 789/3f/3 5 48721/? 2 /? 3 /3 4 



56a/7 8624a/VL 
1658135/3 2 /3 4 /3 5 



3/3g 19557^/35 



2186184V55 
, 27147/3 3 /3| 



101640V35 

t 35403/? 2 /3 5 /3 6 
11440V715 
657095/3 3 AtAs 



88\/7 104104VTT 220220\/7 528528a/9T 



141723/3 3 ^I 



110793a/7/3 3 /3 2 



245625/3f/3 5 



1041040a/7 

| 46892/3 4 /3 5 /3 6 
36465a/143 



+ 



5348200 
327/Sf 



1457456VH 

64461/3 5 /3 2 " 



.-)74(i\ I I ! 369512 v VI 



(29) 



A20 — 



; 3 A 2 / 3 Z 

7T 3 / 2 Q 



297/3 2 /3 3 | 20277/3 2 /3 5 
2240a/7 68992VTT 



80 1 81/3 2 /3 3 /3 4 2047545^ 2 /3 3 /3 6 16455195/3 2 /3 4 /3, 



54208V35 
252207/3 2 /3 5 /3 6 
18304a/715 



+ 



352352a 



+ 



5829824V55 



81/3| 56025/3|^ 5 177669/3 3/ 3f 



704a/7 75712a/vI 352352x77 



8432595/3 3 /3 4 /3 6 1113129/3 3/ 8g 
1409408\/9l 1665664a/7 

9802305/3|/3 5 



1037259y / 7/3 3 /3g 
8557120 
29 9 061/3 4 /3 5 /3 6 



11659648V11 
31455/3f 



38896VV 
3679965/3 5 /3 2 



91936VH 2956096VH 



(30) 



are the volume and surface second-order contributions re- 
lated to (p°(r) contribution (see Eqs. (17), (18), (20), (24)), 
and 



A21 — —- 



5 AV 3 Z 
7rJ 2 ro 



176933/3 4 /3, 



477/3 2 /3 3 , 3719/3 3 /3 4 



15680a/35 194040a/7 
627219/3 5 /3 6 



6342336VH 5725720a 



(31) 
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Fig. 1. Experimental and theoretical values of the PEDM as 
well as macroscopic and microscopic contributions into the 
PEDM for Th isotopes. Points are experimental data from 
Refs. 13,25,26 . Details for theoretical lines see in text. 



A 



nJQro 



459/3 2 /3 3 , 9623/3 3 /3 4 



7840 



32881/3 4 /3 5 



35 258720a/7 
702081/3 5 /3 6 



604032VH 3271840a 



(32) 



are the volume and surface second-order contributions con- 
nected to v? 1 (r) contribution (see Eqs. (17), (18), (19), 
(25)), r = Rop/A 1 / 3 . This expression of the PEDM is 
obtained with the help of symbolic computation software 
Mathematica. 

We propose at evaluation of the PEDM that ratio of 
potentials (p 1 (r) / (p° (r) is the same order as /3g, therefore 
we take into account terms proportional to the product 
of deformations Pi(3t in A21 or A21 and neglect by the 
next order terms /3ef3ei/3e". This proposal is natural for 
the hierarchy of solutions in the form of the perturbation 
series. 

The first term in Eq. (27) was obtained in Refs. [2|3], 
Eqs. (27) and (28) were derived in Refs. [HUH], an d Eqs. 
(29)- (32) are found for the first time. 



3 Discussion and conclusions 

The total value of PEDM Aot is the sum of macroscopic 
-Dmacro and microscopic shell-correction Anicro contribu- 
tions [9"1 H4"1IT2] calculated for the same shapes of the pro- 
ton and neutron surfaces [9], i.e. 



Ant = A 



A 



(33) 



The macroscopic part of PEDM can be evaluated by using 
Eqs. (26)-(32). 

The total values of PEDM evaluated in the frame- 
work of various models are compared with the experimen- 
tal data for thorium isotopes in Fig. 1. The experimen- 
tal data are taken from Refs. [13,25,26]. Our calculation 
of the macroscopic part Anacro is done with the help of 
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Fig. 2. Dependencies of the total macroscopic PEDM eval- 
uated in the first Di maC ro and second -D macro orders on 
quadrupole P2 deformation as well as the same dependencies 
of contributions D v \, D V 2o, D V 2i, D s i, D S 2o, and D S 2i into the 
PEDM. The quadrupole and octupole deformations of 220 Th 
are only taken into account. 



Fig. 3. Dependencies of the total macroscopic PEDM evalu- 
ated in the first Di macro and second -D macr o orders on octupole 
Pz deformation as well as the same dependencies of contribu- 
tions D v i, D V 2o, D V 2i, Dsi, D S 2o, and D S 2i into the PEDM. 
The quadrupole and octupole deformations of 220 Th are only 
taken into account. 



Eqs. (26)-(32) using recent parameter values of the droplet 
model J = 32.5 MeV, Q = 29.4 MeV, r = 1.16 [§]. 
The values of multipole deformation parameters fit and 
microscopic part of PEDM D m i cro are taken from Ref. 
[T4"] . The results obtained in our model well agree with 
the experimental data for 220-228^1^ gee pjg l. The to- 
tal values of PEDM calculated by Butler and Nazarewicz 
using the droplet model approach for macroscopic part 
Act bn = dmacro dm + -Dmicro .13] are shown in Fig. 1 
too. The results of the droplet model approach under- 
estimate the experimental data for 223 - 228 Th, see Fig. 
1. The values of PEDM obtained in the framework of 
cranking Skyrme-Hartree-Fock approach Dshf [E] are 
also presented in Fig. 1. The values of PEDM evaluated in 
the cranking Skyrme-Hartree-Fock model underestimate 
the experimental data for 222 — 227 Th anc i overestimate the 
ones for 229 . 230 Th. Comparison of the PEDM values calcu- 
lated in the framework of various models with the exper- 
imental data for thorium isotopes in Fig. 1 suggest that 



our proposal, that the proton and neutron surfaces are 
geometrically similar (see Eq. (2)), is reasonable. 

The values of macroscopic part of PEDM evaluated 
in our model -D ma cro are larger than the ones obtained in 
the framework of the droplet model D macro DM, see Fig. 
1 and Refs. [6j[7l[8j[9] . Comparing the results in Fig. 1 we 
conclude that the microscopic contribution to the PEDM 
-Dmicro is smaller than the macroscopic one -D macro . 

Let us study the role of second-order contributions into 
the macroscopic part of the PEDM in well-deformed nu- 
clei. We consider nuclei with quadrupole and octupole de- 
formations at the beginning. Here we neglect by the mi- 
croscopic part of PEDM for simplicity. Dependence of the 
PEDM on the quadrupole deformation value at fixed value 
of octupole deformation is presented in Fig. 2, while de- 
pendence of the PEDM on the octupole deformation value 
at fixed value of quadrupole deformation is presented in 
Fig. 3. Fixed values of quadrupole or octupole deforma- 
tions are pointed in Figs. 2-3. These values of the defor- 
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—- o— 



0,0 0,2 0,4 P 2 0,6 

Fig. 4. The same as in Fig. 2, but the high-order multipole de- 
formations Pa, Ps and p 6 of 220 Th are also taken into account. 




1 1 i 1 -r- 

0,0 0,2 0,4 p 0,6 

Fig. 5. The same as in Fig. 3, but the high-order multipole de- 
formations Pa, Ps and /3 6 of 220 Th are also taken into account. 



mation parameters are typical for the ground-state of re- 
flection asymmetric actinides |14) . 

The macroscopic PEDM consists of six contributions 
Art) A20, A21, Ai, A20, and D s2 i, see Eqs. (26)-(32). 
Contributions of all these terms into the PEDM as well 
as the total first and second orders macroscopic PEDM 
values are shown in Figs. 2-3. Comparing various lines in 
Figs. 2-3 we conclude. 

- The total first order contribution of the PEDM -Di maC ro 
D v i + D s i is mainly determined the value of PEDM at 
small values of the deformation parameters. The influ- 
ence of the second order terms rises with the values of 
the deformation parameters. 

- Contributions At, Ai, A 1 macro, A21, and D s2 i of 
the PEDM depend on variable deformation parame- 
ter linearly, while Anacro, A20 and A20 depend on 
variable deformation parameter quadratically. 

- Surface contribution of any type is approximately twice 
smaller than the volume contribution of the same type. 
(Note that this conclusion depends on the ratio be- 
tween values J and Q. We would remind that we use 
the same values of J and Q as in the droplet model 
at PEDM evaluation, but some authors [f51 [T4"l[l"2] use 



other values J and Q at description of the PEDM ex- 
perimental data.) 

- Terms D v2 \ and D s i\ related to the Coulomb poten- 
tial correction ip 1 give negative contributions into the 
PEDM, while any other contributions are positive. 

- The absolute values of terms A21 and D s2 \ are similar 
to the ones for terms D v2 q and D s2 q. This is support 
our proposal on the hierarchy of solutions of the per- 
turbation series. 

- The contribution of second order terms A20 + A21 + 
A20 + A21 gives small correction for the case of vari- 
able quadrupole and fixed octupole deformations (see 
Fig. 2). However this contribution noticeably enhances 
the value of PEDM at large octupole deformation in 
the case of variable octupole and fixed quadrupole de- 
formations (see Fig. 3). The total PEDM evaluated 
at large octupole and fixed quadrupole deformations 
is large than the one for large quadrupole and fixed 
octupole deformations. 

Nuclei with reflection asymmetry have also non-zero 
values of high-order multipole deformations [T2lfT4] . De- 
pendencies of macroscopic PEDM on the value quadrupole 
and octupole deformation at fixed values of other defor- 
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mations are presented in Figs. 4 and 5 respectively. Val- 
ues of fixed high-order multipolc deformations parameters 
pointed in Figs. 4-5 are typical for the ground-state of re- 
flection asymmetric actinides [14] . 

The values of high-order multipole deformations are 
smaller than the ones for quadrupole or octupole deforma- 
tions as a rule [12lll4j . Nevertheless high-order multipole 
deformations enhance the value of PEDM noticeably, com- 
pare results presented in Figs. 2-5. The influence of second 
order terms is also strengthening by high-order multipole 
deformations. 

Qualitatively similar results related to the influence 
of second order contribution on the PEDM was obtained 
numerically in Ref. [12] , but contributions D V 2i and D S 2i 
related to tp 1 were skipped in this work. Note that Skalski 
[T2"] evaluate the macroscopic part of the PEDM for the 
case of constant neutron skin (the droplet model approach) . 
As pointed in the introduction the PEDM evaluated in 
the framework of droplet model has additional contribu- 
tion related to the difference between the center of mass of 
neutron skin of uniform thickness and the center of mass 
of the nucleus, which reduces the value of the PEDM in- 
duced by density redistribution [B] , compare also lines de- 
noted as -Dmacro and Dmacro dm in Fig. 1. This neutron 
skin contribution into the PEDM is zero in the case of 
neutron skin formed by geometrically similar neutron and 
proton surfaces. Therefore numerical comparison of our 
and Skalski's results cannot show the difference between 
the first and second order contributions into the PEDM 
obtained in various approaches, because the comparison 
shows mainly the difference related to the neutron-skin 
contribution. 

We evaluate the PEDM of the hyperdeformed state 
of 152 Dy. The values of deformation parameters of 152 Dy 
in hyperdeformed state are (3 2 = 0.61, (3^ — 0.1, (3 4 = 
0.11 P 5 = 0.05 and fa = [12] . The values of macro- 
scopic part of PEDM obtained in the first and second 
orders using Eqs. (26)-(32) are D ma cro = 0.67 e fm and 
-Dimacro = 0.46 e fm correspondingly. The microscopic 
shell-correction part of PEDM evaluated for geometrically 
similar proton and neutron surfaces is -D s hcii = —0.34 e 
fm [T2J . As the result the total values of PEDM found in 
the cases of applying the first and second orders calcu- 
lation of the macroscopic part of PEDM are -D maC ro + 
Ashen = 0.33 e fm and Di maC ro + D s hd\ = 0.12 e fm 
respectively. Note that the total values of PEDM eval- 
uated using the exact numerical calculation of the macro- 
scopic contribution in the framework of the droplet model 
is £> macro skalski + -Dshdi = 0.06 e fm [T2J. Hereby, the 
PEDM depends strongly on the second-order terms in 
well-deformed nuclei as well as on the neutron skin shape. 

In conclusion, the expression for macroscopic PEDM 
taken into account the first and second order terms on the 
parameters of multipole deformations is obtained in the 
case of geometrically similar proton and neutron surfaces 
of reflection asymmetric nuclei. 

The second order terms are important at large values 
of deformation parameters. The second order terms are es- 



pecially important for nuclear shape with non-zero values 
of high-order multipole deformation parameters. 

Second order contributions D V 2o and D S 2o enlarge the 
value of PEDM obtained in the first order. In contrast 
to this contributions D V 2i and D S 2i decrease the value of 
PEDM obtained in the first order. Compensation of these 
contributions into the PEDM value occurs at small values 
of deformations. However contribution D V 2q + -D S 20 into 
the PEDM value is larger than contributions D V 2i + -D S 2i 
at large values of deformation parameters. 

The reduced probabilities of dipole transitions B(E\) 
are measured in various experiments 13,22,23,24, 251121)] . 
Note that B(E1) is proportional to the squared value of 
the PEDM D 2 [Tl[T71[IHl[Tl[2n][2J]. Therefore the second 
order contribution into PEDM leads to significant varia- 
tion of B(E1) in well-deformed nuclei. 

The obtained expression can be easily applied for esti- 
mation of the macroscopic PEDM and strength of dipole 
transition probabilities in various well-deformed nuclei. 
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